function res=thetaR(t)
    res=0.01*t;
end